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SUMMARY 


Nonlinear equations of motion are derived for the Rim Inertial 
Measuring System (RIMS), which- is mounted in a strapped-down configura- 
tion on a carrier vehicle. The RIMS can be used for measuring angular 
rates and linear accelerations and the equations derived can be easily 
interfaced with the dynamic model of a user-defined carrier vehicle. 

Two methods are presented for rate measurement, and the results of the 
nonlinear simulation are presented. 

INTRODUCTION 

The Rim Inertial Measuring System (RIMS) consists of very small 
annular momentum control devices (AMCD's), mounted in a strapped down 
configuration, which are used to measure angular rates and linear accel- 
erations of a moving carrier vehicle such as an aircraft or a space- 
craft. The concept of RIMS was introduced in reference 1, and that of 
AMCD was introduced in reference 2. An AMCD consists of a spinning rim 
suspended in a set of electromagnetic actuators (suspension stations) . 

The rim is spun by a noncontacting electromagnetic spin motor. As de- 
scribed in reference 1, a minimum of two AMCD's will be required for 
measuring angular velocities about three orthogonal axes. However, a 
RIMS consisting of a single AMCD would be sufficient for the purpose of 
analyzing the RIMS performance. 

In this report, nonlinear equations of motion are derived for the RIMS. 
The equations can be easily interfaced with a carrier model and need 
only the following user-supplied information to be complete: (1) nominal 

location of RIMS center relative to the carrier center of mass (c.m.), 

(2) carrier angular velocities and linear velocities in the carrier 
coordinates, and (3) carrier attitude and position in order to initialize 
the equations. The carrier angular and linear velocities must be obtained 
from the user's vehicle model. In this report, three actuator stations 
(fig. 1) are assumed. Each actuator station consists of an axial actuator 
and a radial actuator. It is also assumed that actuator forces do not 
affect the vehicle motion. Although infinite bandwidth actuators are 



assumed in the derivations, actuator models can be easily inserted for a 
more realistic simulation. Results of the computer simulation of the 
RIMS are presented. 

SYMBOLS 


A, B B , B 
’ f n^ n 2 

C 


C ( ) 
P 


f f 
r . z . 

J J 

f ,f 
r . z . 
cj cj 

f 

vz 


g 

®m 


g 


m 


H 

H 


coefficient matrices 

measurement matrix 

cross-product matrix of a vector 

unit vector in k-direction (k = 1,2,3) 

actuator force vector at station j 

radial and axial actuator forces at station j 

radial and axial actuator command forces at station j 

z-direction force acting on vehicle 

axial actuator force vector 

feedback gain matrix 

gravitational acceleration vector 

measured axial centering errors ( 3 x 1 vector) 

measured axial centering error rates (3x1 vector) 

z-axis gravitational acceleration 

rim angular momentum about spin axis 

rim angular momentum vector 

transverse-axis rim inertia 

rim inertia matrix 

vehicle inertia matrix 

Kalman gain matrix 
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k d ,k d 

a r 
K p ,Kd 
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r 


m 


v 


n 


v 


r 



v 

a 


v 

r 


w 


g 

(X.Y.Z) 


X 


Xl 


(x , Y ) 
v r r 
v v 

z 

r 


rate gains 

proportional and rate gain matrices 
proportional gains 
moment acting on rim 

kth component of moment acting on rim of unit A 
kth component of moment acting on rim of unit B 
rim mass 
vehicle mass 

white noise torque acting on vehicle 
rim radius 

coefficient matrix defined in equation (29) 

matrix consisting of the first two columns of T 

transformation matrix (j = 1,2,3) 

spin motor torque 

torque acting on vehicle 

jth column of T 

covariance intensity matrix of n v 
actuator noise 

inertial velocity of rim center of mass 
proximity sensor noise 

denotes coordinate frames; subscripts are described 
the text 

state vector 

vector defined in equation (52) 

X and Y components of . p 

rim center-of-mass position in the inertial system 



a >“u 

a’ b 

5 ,6 
a. ’ r. 

J J 

6 ,6 

a . r . 
cj cj 

6 ,6 
a . r 
rcj rcj 


a. 

J 


P j 


4 ,<t> 

e l e 2 


e ,e 
a i a 2 


0 ,0 ,0 
vi v 2 v 3 


ft 

ft 

0) 

c 

0) 


0) 


vehicle center-of-mass position in the inertial system 
coefficient matrices 

axial and radial centering errors at station j 
axial and radial centering error commands at station j 


axial and radial centering error rate commands at 
station j 

defined in equation (36) 

position of the point on the rim closest to actuator 
station j (in a-system) 

position of actuator station j in v-system 
position of rim center-of-mass in inertial system 
position of nominal rim center location in v-system 
position of vehicle center-of-mass in inertial system 
spin motor torque vector 

rotations (in order 1-2) for obtaining e-system 
from v-system 

rotations (in order 1-2) for obtaining a-system 
from inertial system 

vector defined in equation (33) 

= (e v ,e ) T 

vi v 2 

vehicle attitude angles 

rim rotation angle relative to a-system 

matrix defined in equation (35) 




rim spin speed relative to v-system 
absolute angular velocity of a-system 
components of u> (j = 1,2,3) 
absolute angular velocity of rim 
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0) 

r . 
J 

components of oj 


0) 

sp 

= (0,0, a> o ) 


0) 

spin 

= (0,0, * r ) 


(i) 

V 

absolute angular velocity of v-system 

w 

0)’ 

V 

absolute angular velocity of v-system, expressed in 
e-system 


C 

= (Y ,-X ) 

rv rv' 


Superscripts 



T 

transpose of a matrix 


• 

first derivative 


•• 

second derivative 



estimated value of a variable 


MATHEMATICAL MODEL DEVELOPMENT 
Coordinate Systems 

Let (X., Y. , Z.) be the inertial coordinate frame and let (X , Y , Z ) 
111 v a a a J 

denote a coordinate frame obtained by a rotation 0 about the X.-axis 

a l i 

and 0 about the new Y.-axis. The origin of the a-system is the rim 

center of mass which is also the rim center. Let (X ,, Y , , Z , ) denote 

v rb’ rb’ rb' 

a coordinate system fixed to the rim. This system is obtained from the 

(X , Y , Z ) system by a rotation rp (rim roll angle) about the a-axis. 

Let (X y , Y^, Z ) denote a coordinate system which is fixed to the vehicle 

and is centered at the vehicle center of mass. The v-system is obtained 

from the i-system by rotations 0,0,0 in the order 1-2-3. 

Vi’ v 2 ' V3 

-i Let (X^, Y^, Z^) denote a vehicle-fixed coordinate system parallel to the 

v-system, centered at the nominal location of the rim c.m. 

* 

Equations of Motion of the Rim 

The equations of motion are similar to those derived in reference 3 
for the AMCD. 
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Rotational equations. - Let denote the absolute angular velocity 

a. 

of the a-system, and w the angular velocity of the rim with respect 

to the a-system. Then the absolute angular velocity of the rim expressed 
in the a-system is given by 


It} = 0) + It) 

r a spin 


CD 


where 


(0 . = ( 0 , 0 , ) 

spin v 


( 2 ) 


Rim angular momentum is given by 


H = Iu) = I(a) + w .) 
r r r r a spin 


H + to x H = M 
r a r r 


(3) 

(4) 


where is the total momentum acting on the rim. Therefore, 


(d = - 1 


-1 


C(oj)Io) - V* p xf. 
parr £-> a. j 

j = l 3 


- T 


(5) 


where C ( ) denotes the cross-product matrix of a vector, p denotes 
p a. 

the position of the point on the rim closest to actuator j, f. is the 

force at actuator j in the a-system and denotes the spin motor 

torque vector. 

Since the a-system is obtained by rotations 0 and 0 about 

a! a 2 

the X. and Y axes, 
l a ’ 


<d 


a l 


a 2 


L a 3 J 


cos 0 0 

a 2 


sin 0 


a l 


a 2 


L 0 a 3 . 


( 6 ) 
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Therefore, 


4 


1- -1 


- 

“ 

• 

0 

a l 


sec 0 

a 2 

0 

• 

0 

a 2_ 


0 

1 _ 


01 

a l 


l VI 


Since 


to =oi +oi. =oi +ii 
r 3 a 3 spin a 3 r r 


(7) 


( 8 ) 


substituting for oi from equation (6) gives 

cl 


i = -sin 0 0 +o) 

r a 2 a x r 3 


(9) 


Assuming that the magnetic actuator forces act only axially and radially, 
there is no torque generated by the actuator forces about the Z & -axis. 
Assuming further that the spin motor generates a torque T s only about 
the Z -axis, t in equation (5) is given by 

3. S 


T s = (0, 0. T s ) T 


( 10 ) 


(The spin-motor torque is included here for completeness, and is not used 
in the simulation described later.) Thus, oi is obtained by solving 
the differential equation (5) ; the rim spin speed is then given 

by equation (9) . 

Translation equations . - Let v^ denote the inertial velocity of the 
rim center of mass expressed in the a-system. Then 

m (v + C (oi ) v ) 
r \ r p^ a r/ 

where m^ is the rim mass and m^g is the gravitational force (in 
a-system) . 


•t 

j = l 


f . 
J 


m r g 


(ID 
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Let denote the position of the rim c.m. in the i-system. Then 

K • T ? T (%) T 2 T ( 0 a 2 ) v r < 12 > 

where denotes the transformation matrix corresponding to rotation 

about axis k (k = X,Y,Z). 

Expressions for the Actuator Gaps 

The vehicle on which the RIMS would be installed could be a spacecraft, 
aircraft, etc. Let denote the position of the vehicle c.m. in the 

i-system, and w its absolute angular velocity in the v-system. 

It is necessary to obtain expressions for the magnetic actuator gaps. 

An actuator gap or centering error is defined as the distance between the 
point midway between the actuator pole faces and the point on the rim 
(assumed to be a perfect circle of zero thickness for gap calculation) 
which is closest to it. Thus the exact expression for an actuator gap 
will involve (1) deriving an expression for the rim plane, (2) obtaining 
the point on the plane through which the vector normal to the plane 

passes, such that this normal vector also passes through the point S 
located midway between the actuator pole faces, and (3) obtaining the 
point Q on the rim which lies on the line joining the rim center to 
P. The distance between the points S and Q then gives the actuator 
gap. Although the expression obtained in this manner will be exact, it 
is also quite cumbersome to evaluate. However, since the rim motion 
within the actuator pole faces is small, it is advisable to obtain an 
approximate expression for the actuator gaps. An approximate expression 
(in inertial coordinates) for the gap at actuator j is 

S i * P r - \ - T ‘ T ( 0 J T * T ( 0 J T 3 T ( e v 3 k 

♦KkV K)-V( % ) T 2 T (ej[ T 3 ’(e Tj ) p . 
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where p_. denotes the location of actuator j in the v' -system, and 

P r denotes the nominal position of the rim c.m. in the v-system. 
v 

In the vehicle coordinates, the gap at actuator j is given by 


' V - *s(\> (c V2 ) t> (e Vi ) 5j (14) 

The axial centering error at actuator j is 

5 a = e 3 T 6 V ’ j = lj 2 ’ 3 (15) 

3 j 

where e^ denotes a 3 x 1 vector with all zeros except a "1" in the kth 
position. The radial centering errors are obtained by appropriate rota- 
tions about the Z-axis : 


ef T 3 (60°) 6 


vi 


*2 


ej T 3 (-60°) 6 


r 3 


T . 

-ej 6 


v 3 


Radial centering errors are defined as being positive outwards and 


(16) 

(17) 

(18) 


T 3 (<» = 


COS 

-sin 4> 
0 


sin 4 > 0 

cos <j> 0 

0 1 


The position of actuator j 


in the a- system 



(19) 


( 20 ) 
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Let f and f (i = 1, 2, 3) denote the radial and axial actuator 
r i z i 

forces generated at station i. Forces at station i (in the 

v-system) are given by 


f 

r l 

cos 

(60°) 


•*> 

f 

r l 

sin 

(60°) 

(21) 



f 

Zl 




f 

*2 

cos 

/ — \ 

ON 

o 

o 

V / 

1 



CM 

4^ 

sin 

(60°) 

(22) 



f 

Z2 






(23) 



Since the forces f^ are in the v-system, they must be transformed into 
the a-system before their use in equation (5) . 


The control laws for the magnetic actuators are assumed to be pro- 
portional plus derivative type. The radial and axial command forces at 
actuator i are given by 



(24) 


( 25 ) 


10 



where 6 and 6 denote the radial and axial centering error 
r • & • 

Cl Cl 

commands, and 6 , 6 denote the centering error rate commands 

rci rci 

at the ith actuator. These commands are normally zero. K , K_ 

P j °j 

(j = a, r for axial and raidal) are the proportional and derivative 
gains. 

ANGULAR RATE MEASUREMENT SCHEMES 

In order to gain insight into the angular rate measurements, it is 
convenient to examine the linearized equations of motion of the carrier 
vehicle/RIMS. The linearization is performed about nonrotating carrier 
axes. 


Linearized Equations of Motion 
The linearized equations of motion of the rim are given by: 


I 

0 + H 

• 

0 

a 

a l 

a 2 


• • 

• 

I 

0 - H 

0 

a 

a 2 

a l 



(26) 

(27) 


m z = t 3 ^ f + m g 
r r 3 z- r & z 


(28) 


where I is the transverse rim inertia, z the Z-axis rim c.m. 
a r 

position, m g the Z-axis gravitational force, f the 3x1 axial 

T Z Z 

force vector, H is the rim angular momentum (about spin axis) and 
t. is the ith column of matrix T: 

l 



/3r/2 

-r/2 

1 


T = 

- /5r/2 

-r/2 

1 

(29) 


0 

r 

1 



where r represents the rim radius. 
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Linearized two-axis vehicle equations are 


= I Yt + n ) 

V \ V V/ 


( 30 ) 


v 2 


where I and denote the vehicle inertia matrix and the torque 

acting on the vehicle and n v is a disturbance input torque which is 
assumed to be a zero-mean white noise with covariance intensity 
matrix V. It is assumed that RIMS actuator forces have a negligible 
effect on the vehicle. The vehicle c.m. position z^ (Z-axis) is given by 


m z = f + m g 
v v v v b z 
z 


(31) 


where f is the (Z-axis) force acting on the vehicle. Let X r and 
v z v 

Y denote the X and Y components of p t the nominal position of the 
r v r v T 

rim c.m. in the v-system, and let ? = (Y^ , -X^, ) . 

v v 

Dividing equations (26) and (27) by I and subtracting equation 

3 . 

(30) from the result, we have 


M IT 1 

'0 +Q0*+f2G = •=— T f - I _i 

e e v I a z v 
a 


( T V * \) 


(32) 


where 0 denotes 
v 


(0,0V 
V Vi’ v 2 / 


and 


0 = 
e 


0 


0 

- 0 

e l 


a l 

V 1 

0 


0 

- 0 

L e 2 J 


L a 2 

V 2 J 


( 33 ) 
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Jr 


T a - [tl t 2 ] 



H/I a 

0 


(34) 


(35) 


Let 

e = z r " ( z v + C36) 

Dividing equation (28) by m a> equation (31) by m v , subtracting the 
latter from the former and using equations (36) and (30) yield ; 



(37) 


The measured axial centering errors at the three actuator stations 
are given by: 


m 


T 0 + t 3 

a e 3 


e + w 


(38) 


where w^ represents the position sensor noise. The axial actuator 
forces are given by: 


f 

z 


K g 
p & m 


K D g m 


(39) 


where and Kp denote the 3x3 position and rate gain matrices 

and g denotes the measured (or derived) axial actuator centering 
nix* 

error rate vector. 


Estimation Schemes 


Scheme I - direct inversion . - The actuator gains and Kp are 

selected in such a manner that the centering errors are quickly corrected 
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and returned to zero in steady state. The closed-loop bandwidth of the 

centering error loop is usually much higher than that of the vehicle 

motion: that is, the actuator centering errors quickly attain steady 

** • 

state. Thus, substituting 0 = 0 = 0 in equation (32) and solving 

6 6 

for 0 y yield an estimate of the two-axis vehicle angular rates: 


0 =2= Si" 1 

V 





(40) 


The second term in the brackets denotes an error term which is a charac- 
teristic of rate gyros. This error term will be assumed negligible in 
this estimation scheme. The first term can be simplified by using 
equation (39) and making the centering error rate g m equal to 0. 

Thus an estimate of vehicle angular rates based on position sensor 
outputs is given by 


0 = f- 

v I 

a 


T 

T K g 
a p to m 


(41a) 


The same analysis was presented in reference 1 using the frequency- 
domain approach. From equations (31) and (37), the linear acceleration 
(Z-direction) is given by 


1 T- ^ 

Z — 1 3 f + g 

v m 0 z b z 
r 


or 



1 ^ 

— t 3 K g 
m 0 pm 
r r 


(41b) 


Although these rate estimates are based on the linear analysis, it 
is shown below that they are also valid for the general nonlinear equations 
when the vehicle is rotating about all three axes. 
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Nonlinear analysis. - Let (X., Y_, Z_) denote a coordinate system 

6 6 6 

obtained from the v-system by Euler angle rotations <t> and <p in 

e l e 2 

the order 1, 2. The origin of this coordinate system is fixed to the 

rim center and translates with it. Thus the (X g , Y ) plane defines 

the rim plane. Let u) Q denote the angular velocity of the e-system 

relative to the vehicle, expressed in the e-system. Let u)^ be the 

angular velocity of the rim relative to the e-system. Then to = 

T S P 

(0, 0, co q ) where co q = rim spin speed (relative to the vehicle). The 

absolute angular velocity of the rim, expressed in the e-system, is 


co = a)' + co + 0 ) (42) 

r v e sp 

where co^ = absolute angular velocity of the vehicle expressed in the 
e-system. The rim rotational equation is given by [from eq. (4)]: 


( to 1 +co + co ) + C (co' + co ) I (co' + (o + co \ = M (43) 

\ v e sp/ p\v e/r\v e sp/ r 


In steady state (of the magnetic actuator loops) , the e-system does 
not move relative to the vehicle, i.e., co g = to e = 0. Also, spin rate is 
assumed to be constant with respect to the vehicle. Therefore, in 
steady state. 


I fco' J + C fto’^ I (co* + co ) = M 
r\v/ P \ v/ r \ v sp/ r 


(44) 


The second term on the left-hand side is 


- 



- 


- 


0 

-0) ' 
v 3 

0) 1 
v 2 


I to' 
a Vj 


(0 ' 
v 3 

0 

-O) T 
V 1 


I to' 
a v 2 


-> N 

3 

1 

H 
- > 
3 

0 


21 fco' + co ^ 

L a V v 3 o/J 


(45) 
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where I denotes the transverse-axis rim inertia. Therefore, in steady 
state. 


I fu)' + 203 J oj' 
a \ v 3 0 / v 2 

-I /a) 1 +" 2(o ^ a) 1 

a \ v 3 o) v x 

0 


M - 
r 


I to* 
r v 


(46) 


The’ acceleration term I a' 

r v 

in all rate gyros. Ignoring this 


represents an error term which is present 
error term. 


2 oj a)'. + to', oj' 


M 

O V2 V2 v 3 


r 

= 

a l 

- 2a) (o' - (o' (o' 


M 

0 V 1 V3 Vj 


- r a2- 


(47) 


where the right-hand side represents the steady-state torque exerted by the 

magnetic actuators. The torque can be expressed as a moment-arm matrix 

times axial and radial force vector, which, in steady state, can be 

expressed as a proportional gain matrix times the axial and radial centering 

error vector. Thus M can be expressed as 

a 



(48) 


where a is a 2 x 6 matrix, 
a 

If the vehicle is not rotating about the yaw-axis, (o' =0 and 

V 3 

equation (47) can be solved for and uH , the solution being 

basically the same as equation (41a). However, if ? 0, it is 
necessary to use at least one more unit for complete rate measurement. 
The equations for the second unit (unit B) , which consists of a rim 
suspended in the plane, can be similarly derived: 
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r 




2 U) CO 1 + Cl)' Cl)' 


M 

0 v 3 V 3 Vi 


r. 


. 

bi 

-2oj a)' - a)' o' 


\ 

O V2 V 1 v 2 


b 2_ 


In steady state. 



# 49 ) 


(50) 


Equations (47) and (49) represent four algebraic equations in three 

Unknowns: o' , o' , and o' . From these equations: 

Vi* v 2 * v 3 M 

2w (w' - o)' ) = ( M + M \ 

°V v 3 V,J I a V r bj r a2 / 

or 


Cl) 


v 3 


Vl 


+ 


1 

21 a) 
a o 




(51) 


Substituting for o' from equation (51) into the top equation of the vector 

v 3 

equation (49) , a quadratic equation in o^, is obtained. The two roots of 

the equation can be computed. In order to determine the root which represents 

Oy , the first root is substituted in equation (51), which gives o^ . 

Using these values of o^. and o^. , o^ is computed from the top 

part of equation (47), and also from the bottom part of equation (49). If 

these two estimates of a) ' are equal (within some tolerance), the root 

v 2 

chosen was the correct representation of . If not, the other root must 
represent co^ . Three-axis rate information can be obtained from two orthog- 
onal units in this manner. Through algebraic manipulations, it should be 
possible to further simplify the solution. 


Scheme II. - This scheme is based on an application of the Kalman-Bucy 

• • T 

filter. Returning to the linearized equations and denoting = (0^ © v ) 
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5 x 1 


( 52 ) 


Equations (30), (32) and (37) can be expressed in the form: 

B-f+B (t + n ) + B f + B, 
f z ni \ v v/ n 2 vz i 


X = AX + 


'V v 

f a 


(53) 


where X is the 10 x 1 state vector consisting of Xj and Xj, B^, 

B , B are appropriately dimensioned input matrices, and v is the 
n j n 2 3. 

actuator noise vector. For the purpose of estimator design, T v is 

assumed to be zero, and f , n and v are the white-noise inputs 
’ vz v a r 

acting on the system. The Kalman-Bucy filter for this system is given by: 


X = AX + B f f z + K (g m - Cx) (54) 

where C is obtained from equations (38) and (52) . The actuator force 
vector is synthesized as 

f = GX (55) 

where G is the gain matrix obtained from equations (38) , (39) and (52) . 

The steady-state filter has constant gain matrix K, and is more 

suitable for practical implementation than the nonsteady-state version. 

The filter generates estimates of X, which also include estimates of 

ft . The sensor and actuator noise covariance matrices are needed for the 
v 

design. The vehicle input torque covariance matrix (V) can be used 
as a design parameter for the Kalman-Bucy filter. Further analysis is. 
necessary in order to evaluate the performance of this scheme using the 
nonlinear equations for a vehicle rotating about all three axes. 
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NUMERICAL RESULTS 


In order to complete the analysis of RIMS, a computer program was 

written for its simulation using the nonlinear equations developed 

[eqs. (1) - (25)]. A single AMCD was used in this simulation. The 

rim mass was 453.6 gm (1 lb), the rim radius was 5.715 cm (2.25 in.), 

rim spin speed was 1,000 rpm, and the rim center was nominally located 

at the vehicle c.m. With zero initial conditions, a torque T y about 

the X axis was used to force the vehicle for 1 sec and to obtain 
v 

a rate of 0.286 rad/sec at the end of 1 sec. The system was unforced 

for 1 < t _< 4 sec, and a torque -T^ was applied for 4 < t < 5 sec, 

after which time the torque was again zero. Figures 2 to 5 show the 

actual and measured vehicle rate, the actual and measured vehicle 

attitude angle (0 ), the rate measurement error, and the axial and 

radial centering errors at actuator station 1. The actuator-loop 

damping ratio used was 2 (at zero spin speed) in all the computer 

runs. In figures 2 to 5, the radial and axial actuator control-loop 

bandwidth used was 30 rad/sec. Figures 6 to 8 show the RIMS 

performance for = 60 rad/sec. The measurement error as well as the 

actuator centering error is much smaller. The actual and measured 

rate for to . = 100 rad/sec is shown in figure 9. As the actuator 
cl 

control loops are made tighter, the RIMS performance improves. However, 

since the actuator gaps (centering errors) also get smaller with higher 

to ,, more accurate proximity sensors are required if the implementation 
cl 

uses gap measurements rather than force measurements. If the rim 
momentum is increased, higher to^ has to be used for getting satis- 
factory performance. This too results in small gaps and will need more 
accurate proximity sensors. Therefore, from the viewpoint of implementa- 
tion, it would be better to use force measurements [in eq. (40)]. In 
that case it would be advisable to use higher momentum and higher 
since the actuator forces to be measured would be large, resulting in a 
more sensitive measurement scheme. It may also be advisable to use an 
integrator in the actuator control loop. This would result in relatively 
easy measurement of force (by measuring the integrator output) , and also 
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maintain the actuator centering errors near zero. The effect of the 
latter would be beneficial in maintaining the actuators in the linear 
region. Further investigation is needed in this area in order to 
evaluate the dynamic response of RIMS using integrators. 

In summary, the numerical results presented indicate that, with 
properly selected rim momentum and o) cl , the RIMS can provide highly 
accurate rate measurements. 

CONCLUSIONS 

Nonlinear equations of motion have been derived for the RIMS using 
one AMCD. The equations can be easily interfaced with a model of a 
carrier vehicle. The complexity of the equations can be easily increased 
by making additions such as the magnetic actuator models. Two methods 
have been presented for the measurement of vehicle angular velocities. 

The first method is obtained from nonlinear equations, while the 
second method uses linearized equations and the Kalman-Bucy filter. 
Results of the nonlinear simulation using the first method indicate 
that RIMS can offer an attractive and accurate inertial measurement 
system. Further analyses will be required in order to completely 
evaluate the second measurement method, and also to evaluate the per- 
formance of both methods in a stochastic environment. 
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Figure 1. RIMS actuator station locations. 
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